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Abstract 

The reaction of collective oscillations excited in the 
interaction between aperiodically growing Jeans-type 
gravity perturbations and stars of a rapidly rotating 
disk of flat galaxies is considered. An equation is de- 
rived which describes the change in the main body of 
equilibrium distribution of stars in the framework of 
the nonresonant weakly nonlinear theory. Certain ap- 
plications of the theory to the problem of relaxation 
of the Milky Way at radii where two-body relaxation 
is not effective are explored. The theory, as applied 
to the Solar neighborhood, accounts for observed fea- 
tures, such as the shape for the velocity ellipsoid of 
stars and the increase in star random velocities with 
age. 

Subject headings: galaxies: kinematics and dynamics — 
galaxies: structure — instabilities 

1. Introduction 

Relaxation of stellar distribution in galaxies is 
not completely understood yet. Lynden-Bell (1967) 
and later Shu (1978) proposed violent relaxation in 
spherical-like (that is, nonrotating) protogalaxies in 
not virial equilibrium. The associated relaxation for 
an individual to gain or lose energy occurs on the ex- 
ceedingly short time much smaller than one typical 
radial period and well before the rotating galaxy disk 
is formed. The relaxation, however, does not stop at 
this stage. There are numerous observations clearly 
showing that there exists ongoing slow relaxation (on 
the time scale of 10 — 20 rotation periods or even 
larger) in the rapidly rotating disk of Milky Way's 
Galaxy (Wielen 1977; Binney & Tremaine 1987, p. 
470; Gilmore, King, & van der Kruit 1990). This 
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slow relaxation of the distribution of young stars 
which were born in the equilibrium disk of the Galaxy 
results in a randomization of the velocity distribu- 
tion ( "Maxwellianization" ) and a monotonic increase 
of the stellar random velocity dispersion ("heating") 
from about 15 kms" 1 for the youngest stars to about 
40 km s _1 for the oldest disk stars with increasing stel- 
lar ages from ~ 10 6 yr to - 4 10 9 yr. Widen (1977) 
has found that the observed increase of the velocity 
dispersion of disk stars with increasing age indicates 
strongly a significant irregular gravitational field in 
the galactic disk. The irregular field causes a rapid 
diffusion of stellar orbits in velocity (and positional) 
space. The nature of this relaxation should be quite 
different from the violent relaxation. Various mech- 
anisms for the slow relaxation have been proposed. 
See, e.g., Grivnev & Fridman (1990) for a review of 
the problem. In the present work, we elaborate upon 
the idea of the collective relaxation: unstable gravity 
perturbations in the disk affect the averaged veloc- 
ity distribution function. The instabilities and subse- 
quent collective relaxation occur near the equilibrium 
and the perturbations remain relatively small which 
makes this process very different from what occurs 
during violent relaxation. 

Apparently, Toomre (1964, p. 1237), Goldreich & 
Lynden-Bell (1965), Barbanis & Woltjer (1967), Kul- 
srud (1972), and Jenkins & Binney (1990) have first 
suggested instabilities as a cause of enhanced relax- 
ation in disk-shaped rapidly rotating galaxies. It was 
stated that because of its long-range Newtonian forces 
a self-gravitating medium (a stellar "gas," say) would 
possess collective properties: collective, or coopera- 
tive motions in which all the particles of the system 
participate. These properties would be manifested 
in the behavior of small gravity perturbations arising 
against the equilibrium background. Collective pro- 
cesses are analogous to two-body collisions, except 
that one particle collides with many which are col- 
lected together by some coherent process such as a 
wave or an unstable perturbation. The collective pro- 
cesses are random, and usually much stronger than 
the ordinary two-body collisions and leads to a ran- 
dom walk of the particles that rapidly takes the com- 
plete system toward thermal equilibrium. 

We present (for the first time as far as we are 
aware) a quantitative theory of relatively slow relax- 
ation on the Hubble time ~ 10 10 yr of self-gravitating, 
rapidly rotating stellar disks of flat galaxies toward a 
thermal quasi-steady state by collective effects. In 



the process a star "collides" with inhomogeneities of 
a galactic gravitational field which result from the de- 
velopment of the Jeans instability. ^ We find that it 
successfully accounts for several basic observations of 
the Milky Way. 

2. Equilibrium 

In the rotating frame of a disk galaxy, the collision- 
less motion of an ensemble of identical stars in the 
plane of the system can be described by the Boltz- 
mann equation for the distribution function f(r,v,t) 
without the integral of collisions (Lin, Yuan, & Shu 
1969): 
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where r, ip, z are the galactocentric cylindrical coor- 
dinates, the total azimuthal velocity of the stars was 
represented as a sum of the random v v and the ba- 
sic rotation velocity V ro t — r£l, v r is the velocity in 
the radial direction, and the epicyclic frequency k(t) 
is defined by k = 2fi[l + {r / 2Q) (dQ /dr)} 1 ^ 2 . The 
quantity f2(r) denotes the angular velocity of galac- 
tic rotation at the distance r from the center, and k 
varies from 2Q for a rigid rotation to £1 for a Keplerian 
one. Random velocities are small compared with rQ. 
Collisions are neglected here because the collision fre- 
quency is much smaller than the cyclic frequency f2. 
In the kinetic equation (1), $(r, t) is the total gravita- 
tional potential determined self-consistently from the 
Poisson equation V 2 $ = 4irG J fdv = AirGn, where 
n is the volume density. 

The equilibrium state is assumed, for simplicity, to 
be an axisymmetric and spatially homogeneous stellar 
disk. The distribution function may also be a func- 
tion of r, for instance, in the case of an inhomogeneous 
disk, in which case the theory is significantly compli- 
cated (Alexandrov, Bogdankevich, & Rukhadze 1984, 



2 The classical Jeans- type instability of small-amplitude gravity 
disturbances is one of the most frequent and most important 
instabilities in the stellar subsystems of galaxies. The Jeans 
instability is driven by a strong nonresonant interaction of the 
gravity fluctuations with the bulk of the particle population, 
and the dynamics of Jeans perturbations can be characterized 
as a fluidlike interaction. Combined with the familiar Lin-Shu- 
Kalnajs dispersion relation this is a venerable suggestion as to 
why fiat galaxies almost always exhibit spiral structure (Binney 
& Tremaine 1987, p. 336). 



p. 425). Secondly, in our simplified model, the per- 
turbation is propagating in the plane of the disk. 
This approximation of an infinitesimally thin disk is 
a valid approximation if one considers perturbations 
with a radial wavelength that is greater than the typ- 
ical disk thickness. We assume here that the stars 
move in the disk plane so that v z — 0. This al- 
lows us to use the two-dimensional distribution func- 
tion / = f(v r ,v tp ,t)5(z) such that J fdv r dv v dz = a, 
where a is the surface density. We expect that the 
waves propagating in the disk plane have the great- 
est influence on the development of structures in the 
disk. The latter suggestion is strongly supported by 
numerical simulations (Hohl 1978). 

The disk in the equlibrium is described by the fol- 
lowing equation: 



dv^ 



= 0, 
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where df c / dt = and the angular velocity of rotation 
fi(r) is such that the necessary centrifugal accelera- 
tion is exactly provided by the central gravitational 
force r£l 2 = d$ c /dr. Equation (2) does not deter- 
mine the equilibrium distribution / e uniquely. For 
the present analysis we choose f c in the form of the 
anisotropic Maxwellian (Schwartzschild) distribution 
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The Schwarzschild distribution function is a function 
of the two epicyclic constants of motion £ = v\j2 and 
fQO(ro), where tq = r + (2fi/ n 2 )v ip . These constants 
of motion are related to the unperturbed star orbits: 

Vj_ 

r = [sin(^o — nt) — sin 0q] ; v r = v± cos (<po — Kt) ; 
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where vx, 4>o are constants of integration, v±/nrQ ~ 
p/rQ -C 1, p is the mean epicycle radius, and we fol- 
low Lin et al. (1969), Shu (1970), and Griv & Pe- 
ter (1996) making use of expressions for the unper- 
turbed epicyclic trajectories of stars in the equilib- 
rium central field Q c (r). In equations (3) and (4), 
ro is the radius of the circular orbit, which is chosen 
so that the constant of areas for this circular orbit 
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r 2 (dipo/dt) is equal to the angular momentum imte- 
gral M z = r 2 (dip/dt), and v\ = v 2 + (2n/n) 2 v 2 . 
Also, ip is the position angle on the circular orbit, 
(dipa/dt) 2 = (l/r )(d® e /dr)o = O 2 . The quanti- 
ties Cl, k, and c r are evaluated at tq. In equation 

(3) the fact is used that as follows from equations 

(4) in a rotating frame the radial velocity dispersion 
c r and the azimuthal velocity dispersion c v are con- 
nected through c r w (2£1/k)c v . In the Solar vicin- 
ity, 20/k ss 1.7. The distribution function f e has 
been normalized according to J_ J_ f e dv r dv v = 
2tt(k/2£1) J °° v±dv±f e = (T , where <r c is the equilib- 
rium surface density. Such a distribution function for 
the unperturbed system is particularly important be- 
cause it provides a fit to observations (Lin et al. 1969; 
Shu 1970). It is this equilibrium that is examined for 
stability. 

3. Collisionless Relaxation 

We proceed by applying the standard procedure 
of the weakly nonlinear (or quasi-linear) approach 
(Galeev & Sagdeev 1983; Alexandrov et al. 1984, 
p. 408; Krall & Trivelpiece 1986, p. 512) and de- 
compose the time dependent distribution function 
/ = /o(v, t) + /i(v, t) and the gravitational poten- 
tial $ = $oM) + with Lh/fol < 1 and 
|$i/$o| "C 1 for all r and t. The functions /i and 
$i are oscillating rapidly in space and time, while 
the functions /o and $o describe the slowly develop- 
ing "background" against which small perturbations 
develop; f Q (t = 0) = f e and $ (* = 0) = $ c . The 
distribution /o continues to distort as long as the dis- 
tribution is unstable. Linearizing equation (1) and 
separating fast and slow varying variables one obtains: 

dfi _ 0*! df , 1 S/o 
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where d/e?t means the total derivative along the star 
orbit (4) and (• • •) denotes the time average over the 
fast oscillations. To emphasize it again, we are con- 
cerned with the growth or decay of small perturba- 
tions from an equilibrium state. 

In the epicyclic approximation, the partial deriva- 
tives in equations (5) and (6) transform as follows 
(Lin et al 1969; Shu 1970; Griv & Peter 1996): 
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To determine oscillation spectra, let us consider the 
stability problem in the lowest WKB approximation: 
the perturbation scale is sufficiently small for the disk 
to be regarded as spatially homogeneous. This is ac- 
curate for short wave perturbations only, but qualita- 
tively correctly even for perturbations with a longer 
wavelength, of the order of the disk radius R. In this 
local WKB approximation in equations (5) and (6), 
assuming the weakly inhomogeneous disk, the pertur- 
bation is selected in the form of a plane wave (in the 
rotating frame): 

fl, *1 - \sf, (5$ {e^rr+imif-i^t + c c ) ^ (g) 

where Sf, (5$ are amplitudes that are constant in 
space and time, m is the nonnegative azimuthal mode 
number, = lu — mCl is the Doppler-shifted wave- 
frequency, and \k r \R > 1 (Griv & Peter 1996). The 
solution in such a form represents a spiral wave with 
m arms whose shape in the plane is determined by 
the relation k r (r — m) = —m{ip — ipo). With ip in- 
creasing in the rotation direction, we have k r > for 
trailing spiral patters, which are the most frequently 
observed among spiral galaxies. A change of the sign 
of k r corresponds to changing the sense of winding 
of the spirals, i.e., leading ones. With m — 0, we 
have the density waves in the form of concentric rings 
that propagate away from the center when k r > or 
toward the center when k. r < 0. 

In equation (5) using the transformation of the 
derivatives d/dv r and d/dv v given by equations (7), 
one obtains the solution 
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where = — oo) — > 0. In this equation making use 
of the time dependence of perturbations in the form 
of equation (8) and the unperturbed trajectories of 
stars given by equations (4) in the exponential factor, 
it is straightforward to show that 
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where J; (x) is the Bessel function of the first kind of 
order I, x — k*v±/n ~ fc*p, fc» = k{l + [(20/k) 2 — 
1] sin 2 iji} 1 / 2 is the effective wavenumber, -0 is the 
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pitch angle of perturbations, tan^ = k^/kr = m/rk r , 
and we used the expansion 

oo 

exp(±i6sin0 o ) = ^ J n (b) exp(±in<f>o) 

n— — oo 

and the usual Bessel function recursion relation 

Jm(x) + ^-i(x) = (2//x)^(x)- 

In the equation above the denominator vanishes when 
— Ik — > 0. This occurs near corotation (I = 0) 
and other resonances (I = ±1, ±2, • • •). The Lindblad 
resonances occur at radii where the field (d/dr)§\ 
resonates with the harmonics I = — 1 (inner reso- 
nance) and I = 1 (outer resonance) of the epicyclic 
(radial) frequency of equilibrium oscillations of stars 
k. Clearly, the location of these resonances depends 
on the rotation curve and the spiral pattern speed; 
the higher the m value, the closer in radius the res- 
onances are located (Lin et al. 1969; Shu 1970). In 
this paper, only the main part of the galactic disk is 
studied which lies sufficiently far from the resonances: 
below in all equations lj* — Ik ^ 0. 

We substitute the solution (10) into equation (6). 
Taking into account that the terms I ^ n vanish for 
axially symmetric functions /o, after averaging over 
4>o we obtain the equation for the slow part of the 
distribution function: 

dfo V- ,,,, ,2 d Kk Z 2 J ; 2 (x) df 
at L - J ov± v±y w» — Ik ov± 

k i = — oo 

(11) 

As usual in the weakly nonlinear theory, in order 
to close the system one must engage an equation for 
$i,k- Averaging over the fast oscillations, we have 

(a/at)|$ lik | 2 = 23^|o> 1 , k | 2 , (12) 

where suffixes k denote the kth Fourier component. 

Equations (11) and (12) form the closed system of 
weakly nonlinear equations for Jeans oscillations of 
the rotating homogeneous disk of stars, and describe 
a diffusion in velocity space. The spectrum of oscil- 
lations and their growth rate are (Griv, Rosenstein, 
Gcdalin, & Eichler 1999a: Griv, Gedalin, Eichler, & 
Yuan 2000a) 

k 2 c 2 r A e - x h{x) 

2vGa \k\ " ^ w.-Ik ' { ' 

i = — oo 



and $suj* t j ~ ^/AnGa^e x I(x)/x <; O, respectively. 
In the Solar vicinity, £1 s=a 3 10~ 8 yr _1 . Here, h{x) 
is a Bessel function of an imaginary argument with 
its argument x ~ kip 2 and p = c r /K is now the 
mean epicyclic radius. A very important feature of 
the instability under consideration is the fact that it 
is aperiodic (the real part of the wavefrequency van- 
ishes in a rotating frame we are using). Usually the 
quasi-linear theory is applied when the growth rate 
is small compared with the real part of the wavefre- 
quency as for the case of the resonance interaction 
w = k • v, where v is the velocity of the particle in- 
volved in the interaction. However, the theory can 
be applied also to aperiodic instabilities (Shapiro & 
Shevchenko 1963; Alcxandrov et al. 1984, p. 420; 
Krall & Trivclpicce 1986, p. 531). A further simpli- 
fication results from restricting the frequency range 
of the waves examined by taking the low frequency 
limit (|aj* | less than the epicyclic frequency of any disk 
star). In the opposite case of the high perturbation 
frequencies, |u>»| > k, the effect of the disk rotation 
is negligible and therefore not relevant to us. This 
is because in this "rotationless" case the star motion 
is approximately rectilinear on the time and length 
scales of interest which are the wave growth/damping 
periods and wavelength, respectively (Alexandrov et 
al 1984, p. 113). Thus, the terms in scries (10)-(13) 
for which \l\ > 2 may be neglected, and consideration 
will be limited to the transparency region between 
the turning points in a disk (between the inner and 
outer Lindblad resonances). In this case, in equations 
(10)-(13) the function A(x) = exp(— x)Ii(x) starts 
from A(0) = 0, reaches a maximum A max < 1 at 
x w 0.5, and then decreases. Hence, the growth rate 
has a maximum at x < 1 (see Griv, Yuan, & Gedalin 
1999b, Fig. 2 in their paper). 

In general, the growth rate of the Jeans instability 
is high ISw*,,]! <~ O; perturbations with wavelength 
A j ~ 2np have the fastest growth rate (Morozov 1981; 
Griv & Peter 1996; Griv et al. 1999a). In the Solar 
vicinity of the Galaxy, Aj = 2 — 4 kpc. 

4. Astronomical Implications 

As an application of the theory we investigate 
the relaxation of low frequency and Jeans-unstable, 
< k and J 2 < 0, respectively, oscillations in 
the homogeneous galactic disk. Indeed, already in 
the 1940s it was observed that in the Solar neigh- 
borhood the random velocity distribution function of 
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stars with an age t ^ 10 8 yr is close to a Schwarzschild 
distribution — a set of Gaussian distributions along 
each coordinate in velocity space, i.e., close to equi- 
librium along each coordinate (Chandrasekhar 1960; 
Ogorodnikov 1965; Binney & Tremaine 1987, p. 471). 
In addition, older stellar populations have a higher 
velocity dispersion than younger ones.0 On the other 
hand, a simple calculation of the relaxation time of 
the local disk of the Milky Way due to pairwise star- 
star encounters brings the standard value ~ 10 14 yr 
(Chandrasekhar 1960), which considerably exceeds 
the lifetime of the universe. According to our ap- 
proach, collisionless relaxation does play a determin- 
ing role in the evolution of stellar populations of the 
Galactic disk. 

Evidently, the unstable Jeans oscillations must in- 
fluence the distribution function of the main, nonres- 
onant part of stars in such a way as to hinder the 
wave excitation, i.e., to increase the velocity disper- 
sion. This is because the Jeans instability, being es- 
sentially a gravitational one, tends to be stabilized 
by random motions (Toomre 1964; Shu 1970; Bertin 
1980; Morozov 1981; Griv & Peter 1996). Therefore, 
along with the growth of the oscillation amplitude, 
random velocities must increase at the expense of cir- 
cular motion, and finally in the disk there can be es- 
tablished a quasi-stationary distribution so that the 
Jeans-unstable perturbations are completely vanish- 
ing and only undamped Jeans-stable waves remain.^] 

In the following, we restrict ourselves to the most 
"dangerous," in the sense of the loss of gravitational 
stability, long- wavelength perturbations, x 2 an d x 2 <C 
1 (see the explanation after eqn. [13]). Then in equa- 
tions (10)-(13) one can use the expansions Jf (x) ~ 
X 2 /4 and e-' x h{x) « (l/2)x - {l/2)x 2 + (5/16)x 3 . 
Equation (11) takes the simple form 

f - 4f • <"» 

where D = (tt/16k 2 ) £ k fe^w* :k | 2 , 3w»„ T > 0, 
and both 9to*,j and <&i,k are functions of t. As is 
seen, the velocity diffusion coefficient for nonresonant 
stars D is independent of v± (to lowest order). This 



3 The age dependence of velocity dispersions for stellar popula- 
tions has always been of particular interest, because the form of 
this relationship allows us to judge whether was any relaxation 
in the galactic disk and even to determine the mechanism that 
was responsible for increase in random star velocities. 

4 In turn, the Jeans-stable perturbations are subject to a reso- 
nant Landau-type instability (Griv et al. 2000a). 



is a qualitative result of the nonresonant character of 
the star's interaction with collective aggregates. 

By introducing the standard definitions dr/dt = 
D(t) and d/dt = (dr/dt)(d/dr), equations (12) and 
(14) are rewritten as follows: 

df a 2 / od 

"5 "ST" = °' "7T - = 23u>*,j, (15) 

OT Ovj_ OT 

which has the solution 

/oK,r) = ^exp(-J). (16) 

(We have taken into account the observations that 
the distribution of newly born stars is close to the 8- 
function distribution, /o( v ±j t = 0) = 5(vj_); Grivnev 
& Fridman 1990.) As is seen from equation (16), 
during the development of the Jeans instability, the 
Schwarzschild distribution of random velocities (a 
Gaussian spread along v r ,v v coordinates in velocity 
space) is established. The energy of the oscillation 
field oc J2k l^i, k| 2 plays the role of a "temperature" 
T in the nonresonant-particle distribution. As the 
perturbation energy increases, the initially monoen- 
ergetic distribution spreads (/o becomes less peaked) , 
and the effective temperature grows with time (a 
Gaussian spread increases): T = 2r oc J D(t)dt oc 

From the above, this mechanism increases the ve- 
locity dispersion of stars in Milky Way's disk after 
they are born. Subsequently, sufficient velocity dis- 
persion prevents the Jeans instability from occuring. 
The "diffusion" of nonresonant stars takes place be- 
cause they gain mechanical (oscillatory) energy as 
the instability develops. The velocity diffusion, how- 
ever, presumably tapers off as Jeans stability is ap- 
proached: the radial velocity dispersion c r becomes 
greater than the critical one c riCr it ~ (2S1/k)ct, where 
ct is the well-known Toomre's critical velocity dis- 
persion to suppress the instability of axial symmetric 
gravity perturbations (Morozov 1981; Griv & Peter 
1996; Griv et al. 1999a). Thus, the true time scale 
for relaxation in the Milky Way may be much shorter 
than its standard value ~ 10 14 yr for the classical 
Chandrasekhar-Ogorodnikov collisional relaxation; it 
may be of the order (Qxt^j) -1 ;> il^ 1 ^ 10 9 yr, i.e., 
comparable with 10 periods of the Milky Way rotation 
in the Solar vicinity. The above relaxation time is in 
fair agreement with both observations (Wielen 1977; 
Knudc, Winther, & Schnedler-Nielsen 1987; Gilmore 
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et al. 1990; Grivncv & Fridman 1990; Meusinger, 
Stecklum, & Rcimann 1991) and TV-body simulations 
of Milky Way's disk (Hohl 1971; Sellwood & Carlberg 
1984; Griv, Gcdalin, Liverts et al. 2000b). 
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